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1 Introduction 

In recent years, the need of accurate theoretical predictions for scattering amphtudes in 
colhder physics requested a strong effort in the development of methods and strategies for 
the calculation of multi-loop Feynman diagrams. In particular, it was recently possible to 
afford many analytic calculations unthinkable up to fifteen years ago. 

This success is largely due to reliable and powerful algorithms, as for instance the so- 
called "Laporta algorithm" [1]. The calculation of a physical observable in perturbation 
theory requires the (numerical or analytical) evaluation of a large number of regularized 
scalar integrals. The Laporta algorithm allows to reduce this large number of scalar integrals 
to a linear combination of a small set of independent scalar integrals, called the "Master 
Integrals" (Mis) of the problem under consideration. The Laporta algorithm is based on the 
Integration-by-Parts Identities (IBPs) [|2J, a set of relations that link scalar integrals with 
a different power of the propagators and of scalar products in the numerator among each 
otheio. The interplay between the "Differential Equations Method" [3] and the techniques 
based on Mellin-Barnes representations of the integrals [5J provides, then, a powerful tool 
for the analytic calculation of the Mis. 

Another important ingredient for the analytic calculation of the higher-order corrections 
to a physical observable is the identification of the base of functions in terms of which the 
Mis can be expressed. This base of functions is strictly related to the structure of the 
thresholds of the Feynman diagrams under consideration. 

The connection between Feynman diagrams with a simple structure of thresholds and the 
functional base of the Harmonic Polylogarithmqj (HPLs) [6j is completely clear. Feynman 
diagrams with a richer structure of thresholds require the introduction of one- fiU\ \TT\ and 
two-dimensional fi2\ [T3| HH |15l |16] extensions of the HPLs, that will be generically referred 
to as Generalized Harmonic Polylogarithms (GHPLs). 

We briefly recall the advantages of using the (G)HPLs, that constitute a well suited 
functional base in which to express the analytic results, i) The structure of the (G)HPLs, 
is strictly connected with the solution, with the Euler method, of the first-order linear 
differential equations satisfied by the Feynman amplitudes, ii) The (G)HPLs constitute a 
base of linearly independent functions. Hi) The base of (G)HPLs provides a perfect control 
on the analytical properties of the Mis and, therefore, of the physical observable that we 
aim to calculate, iv) Finally, there are available numerical routines, that allow a precise 
evaluation of the (G)HPLs in FORTRAN [17], Mathematica [H], C++ [I9]. 

In this paper, we provide a detailed analysis of the GHPLs of a single variable, with 
weights containing square roots. These GHPLs were introduced in [TT] for the analytic 
expression of the Mis concerning the electroweak form factor [201 E]. In [21], ad hoc nu- 
merical routines were made for the evaluation of a subset of functions of this class occurring 
in the calculation of the electroweak NLO corrections to the production of a Higgs boson 
in gluon fusion and its decay in two photons. The purpose of the present analysis is to give 
a general framework for the evaluation of GHPLs with weights containing square roots and 
to show that the evaluation of all of the GHPLs introduced in pjj can be performed using 
already existing numerical routines. 



^Public implementations of the algorithm are available in [3]. 

^Together with the related harmonic [7], nested [8] and binomial sums [9]. 



In section |21 we recall the definition and the basic properties of the 1- dimensional HPLs 
and GHPLs with square roots in the weights. We start on the case of real variable x 
and on the subset of GHPLs that have the following possible weights: {— r, — 4, — 1, 0}. 
Subsequently we discuss their analytic continuation. In section [3l we illustrate how to 
move from the set of GHPLs with square roots in the weights to a set of generalized 
polylogarithms, with linear weights. This is the remark that allows a numerical evaluation 
of the GHPLs with square roots, using the C++ routines of Vollinga and Weinzierl [19] (in the 
following they will be referred to as "VW routines"). In section |U we consider the case in 
which the variable x is complex. We provide a demonstration of the GHPLs scale invariance 
in the complex plain, justifying the use of the VW routines also in this case. In section 
El we introduce additional weights and discuss the transformations of section [3] applied to 
this new extended set. In section [HI we apply the results of this paper to the numerical 
evaluation of the GHPLs involved in the NLO light-fermion electroweak corrections to the 
Higgs boson production in gluon fusion, in the case of complex W and Z masses. Finally, 
in the appendices, we provide the analytic expressions of the linearized GHPLs involved in 
the calculation shown in section O 

2 HPLs and GHPLs of a Real Variable 

In this section, we recall the definition and the properties of the one-dimensional Harmonic 
Polylogarithms (HPLs) of a real variable and their generalization (GHPLs), with square 
roots and linear weights, introduced in [11]. 



2.1 Harmonic Polylogarithms 

The set of functions denominated Harmonic Polylogarithms (HPLs) [B] is defined as re- 
peated integrations of the following three fundamentao "weight functions" : 

f{-l;t) = ^, /(0;t) = i, /(l;t) = ^. (1) 

Note that the functions in Eq. ([1]) have a non-integrable singularity in t = — 1, t = 0, and 
t = 1 respectively. The related HPLs of weight 1 are 

H{-l-x) = j^ ^=log(x + l), (2) 

r df 

H{0;x) = J J =log(x), (3) 

Hil;x) = f-^ =log(l-x), (4) 

where a; is a real variable (a; G M). Since the logarithms have branch cuts on the real axis 
for X < — 1, X < 0, and a; > 1, respectively, the three HPLs in Eqs. (|2|3|^ are real and 
uniquely defined only for x > — 1, x > 0, and x < 1, respectively. Outside these intervals. 



^Note a minus sign in the weight +1 with respect to the Remiddi-Vermaseren definition of [6] 



the logarithms become complex, and a prescription for the approach to the branch cut has 
to be chosen (see section [2. l.ip . 

An HPL with weight 2 or bigger is defined through a repeated integration of the weight 
functions of Eq. ([1]). If w is a vector with w components consisting of a sequence of —1, 0, 
and +1, we define the HPL of weight w + 1 as follows: 

H{a,w,x)= / dt f {a; t) H {w; t) , a = —1,0,1, (5) 

Jo 

with the exception of the case in which the weights are only zeroes, defined as 

i7(0^+i; x)= [ dt /(O; t) H{0^- 1) = [ log^'+i x . (6) 

The singularity structure and analyticity properties of the HPLs derives from the prop- 
erties of the logarithms. A logarithmic singularity in 0, —1 or +1 (and +oo) can occur, 
together with the respective branch cuts in a; > 1, x < 0, and x < —1, as discussed in the 
next section. 

The HPLs satisfy a shuffle algebra according to which a product of two HPLs of weights 
ni and n2 is a combination of HPLs of weight n = ni + n2. Let wi be a vector with rii 
components and W2 a vector with n2 components, both consisting of a sequence of -1, 
and +1. Then we have: 

H{-Wi;x) H{w2;x) = ^ H{w,x), (7) 

w=witiJW2 

where w is a vector with n = rii + n2 components —1, 0, or +1. The sign l±) means that 
the order in the components of Wi and W2 has to be preserved in the sequence w. For 
instance, an HPL with weight 4 can be made out of a product of two HPLs, one of weight 
1 and another of weight 3, or both of weight 2. The relevant formulas are: 

if (a; x) H{b, c, d; x) = H{a, b, c, d; x) + H{b, a, c, d; x) + H{b, c, a, d; x) + H{b, c, d, a; x) , (8) 
H{a, b; x) H{c, d; x) = H{a, b, c, d; x) + H{a, c, b, d; x) + H{a, c, d, b; x) + H{c, a, b, d; x) 

+H{c, a, d, b; x) + H{c, d, a, b; x) . (9) 

The demonstration that HPLs satisfy the shuffle algebra in Eq. ([7]) can be done by induction, 
using integration by parts (see [6]). 

2.1.1 Analytic Continuation 

The HPLs are, in general, complex, depending on the value of the real variable x. In many 
relevant physical cases, the calculation of the Feynman integrals involved in some observable 
is done in a restricted range of x. For instance, if x is related to the squared center of mass 
energy s through the relation x = —s/m?, with m a mass scale of theproblem, the Feynman 
integrals are usually solved in the so-called Euclidean region: a; > Qj. 



*If, moreover, x is related to s through a quadratic relation, as the transformation of variable in Eq. p4l) . 
the Euclidean region is even more restricted: < a; < 1. 



Let us suppose, then, a; > 0. In the range < x < 1 all the HPLs are real. For a; > 1, 
instead, we have a possible cut, corresponding to the HPLs with a +1 in the right-most 
weight, and therefore, ultimately, to the log (1 — x), that has an imaginary part in this 
region. 

In the case of HPLs of weight 1, depending on the prescription adopted, this imaginary 
part is iivr: 

H{l;x) = log(l-x) =\og\l-x\±me{x- 1), (10) 

while H{0; x) and H{—1; x) are real for positive x. 

In the case of HPLs of weight 2, or bigger, an explicit expression for the imaginary part 
can be found, using the shuffle algebra properties to move the weights 1 from the right to 
the left in the sequence and the relation of Eq. f llOp . For instance, for H{0, l;x) = Li2(x), 
the Euler Dilogarithm, we have: 

H(0,l;x) = H{0;x)H{l;x)-H{l,0;x), 

= log(a;) log |1 — x| — if(l,0;a;) ± Z7r6'(x — 1) log(x) . (11) 

The HPL H{1, 0; x) is real for x > 1 and the imaginary part of H{0, 1; x) is explicitly given 
by the last term in Eq. fITT]) . 

Once the searched analytic expression is known in the range of a; > 0, one has to do 
an analytic transformation to move back to the Minkowski region (s > and then x < 0). 
Because of causality, the Mandelstam invariant s has to be assigned a positive vanishing 
imaginary part, s + iO^. Therefore, if x = —s/m?, the case s > is recovered using 

x-^-x'-iO+, (12) 

where now x' = s/iv? > 0. 

If < x' < 1 (— 1 < X < 0), we have to take into account the branch cut connected to 
the weight 0. For HPLs of weight 1 we have 

H{0;x) -^ H{0;-x' -iO+) = H{0;x')-i7i, (13) 

H{l;x) -^ H{l;-x' -iO+) = H{-l;x'), (14) 

H{-l;x) -^ H{-l;-x' -iO+) = H{l;x). (15) 

Therefore, only the log (x) gets the imaginary part. The case of HPLs of weight 2 or bigger 
has to be treated extracting, with the help of the shuffle algebra, the right-most zeroes 
(trailing zeroes). The HPLs that have no zeroes on the right of the sequence of the weights 
do not get imaginary parts. Moreover, moving from x to x', the weights flip in sign. The 
log" X extracted with the algebra are then transformed according to Eq. (II 3p . As a simple 
example, consider the function H{1, —1, 0, x). Using the shuffle algebra we obtain: 

H{1, -1, 0; x) = H{0, -1, 1; x) - H{0, -1; x)H{l; x) + H{1, -1; x)H{0; x) . (16) 

The cut behaviour is explicitly extracted as H{0] x). The other functions, H{1] x), H{1, — 1; x), 
H{0,—l;x), and if(0,— l,l;x) are real for — 1 < x < 0. Therefore, using Eqs. f lT2|) and 



(IT3]) . we have: 

H{1, -1, 0; -x' - iO+) = H{0, -1, 1; -x' - zO+) - H{0, -1; -a;' - iO^)H{l; -x' - iO+) 

+//(!, -1; -x' - iO+)if (0; -x' - iO+) , 
= i/(0, 1, -1; x') - i7(0, 1; x')i/(-l; x') + i7(-l, 1; x')i/(0; x') 
-Z7riJ(-l,l;x'). (17) 

If x' > 1 (x < — 1), we have the superposition of the cuts connected to the weights and 
1. The first cut is taken into account with the analytic continuation discussed just above: 
X — > —x' — iO^ with < x' < 1. Now, if x' > 1, the HPLs with the right-most weight 1 
exhibits an imaginary part, which comes from the log(l — x'). Using the shuffle algebra, 
the weights +1 on the right of the sequence can be moved to the left, and the log (1 — x') 
explicitly extracted. Since s is to be understood with a vanishing positive imaginary part, 
this is also the case of x': x' + iO"*". Therefore, we have: 

log(l-x'-iO+) =log(x'- 1) -ivr, x' > 1 . (18) 

Continuing with the example of if (1, — 1, 0; x), now we start from Eq. fll7p . The cut for 
x' > 1 shows up in the functions -f/^(0, 1; x') and if (— 1, 1; x'). Using again the shuffle algebra 
we can rewrite these functions extracting explicitly the dependence on log (1 — x') as 

iJ(l, -1, 0; -x' - zO+) = if (0, 1, -1; x) + E(\, 0; x')if (-1; a;') - if (1, -1; x')if (0; x') 

-m [if (-1; x')if (1; x') - ff (1, -1; x')] , 
= ff (0, 1, -1; x') + ff (1, 0; x')ff (-1; x') - ff (1, -1; x')ff (0; x') 
-TT^ - m [ff (-1; x) log (x' - 1) - ff (1, -1; x')] , (19) 

Every function in Eq. flT9|) is real for x' > 1 and the imaginary part is explicitly given by 
the last term. 

2.1.2 Transformations of Variables and Numerical Evaluation 

A fast and precise numerical evaluation of the HPLs, for all the values of the real variable x, 
can be done using an appropriate Taylor expansion in the vicinity of a point of analyticity 
of the functions. The strategy is the following (see for instance [T7]): 

1. We focus on the point x = 0. We extract the possible logarithmic behaviour, log"(x), 
of the HPL using the shuffle algebra (we move the rightmost zeroes to the left). The 
HPLs with no zeroes on the right are analytic in x = 0. Correspondingly, each HPL 
takes the form ^„ ^ Pm{x) log"(x), where Pm is a polynomial of degree m. In the case 
X — )• 0~, the imaginary part comes from log'^(x), using the prescription of Eq. (fTOj) . 

2. With an appropriate number of terms in the Taylor expansion, one is able to evaluate 
numerically the HPL in an interval around x = with a given precision. In [17] 
the interval — (v2 — 1) < x < (\/2 — 1) is taken as the central region and using 
Bernoulli numbers and Chebyshev polynomials, the authors evaluate the HPLs in 
double precision using only few terms in the expansion. 



3. Using the properties of the HPLs, one can find suitable transformation formulas for 
the argument in order to map different domains of the real axis back to the central 
value — (\/2 — 1) < x < {\/2 — 1). In so doing, using the formulas found for that 
region, one is able to cover all the possible values of the variable x G R. 

2.2 Generalized Harmonic Poly logarithms 

In some physically relevant cases, it can happen that the weight functions defined in Eq. ([1]) 
are not sufficient to describe the analytic structure of the result. Therefore, additional 
weights (and/or additional structures) have to be taken into account, together with the 
ones introduced in the last section. This gives rise to an enlarged set of functions, called 
Generalized Harmonic Polylogarithms (GHPLs), which maintain the structure and proper- 
ties of the HPLs. 

Let us focus, for the moment, on the GHPLs that are involved in the calculation of 
the NLO light-fermion electroweak corrections to the cross section of production of a Higgs 
boson in gluon fusion and its decay in two photonqfl, as considered in [2T] . This is a restricted 
set, that contains only four weights, denominated as follows: 

G{wi,W2,...,Wn;x), withwj G {-r, -4, -1,0}, (20) 

according to the definitions given below. Let furthermore restrict the analysis to the case 
of real variable x G M (the case of complex x will be treated in section H]). Therefore, we 
consider the following set of weight functions: 

g{-r;x) = . / ^ , (21) 

Va;(x + 4) 

g{w;x) = , with w G {-4, -1,0}. (22) 

X — w 

These functions have an integrable singularity in x = and x = —4, and a non-integrable 
singularity in x = w, respectively. The related GHPLs of weight 1 are 

G{0;x) = log(x), (23) 

G -r; x) = / , = - log \ V , 24 

r dt 

G{w;x) = / = log {x — w) — log {—w) , withw G {— 4, — 1} , (25) 

Jo t-w 

and they have at most a logarithmic singularity in a; = 0, x = —1, —4. G(0; x) and G{w\ x) 
have a branch cut for a: < and x < w, respectively. For these negative values of x, they 
become complex, with imaginary part depending on the prescription of approach to the 
cut. (?(— r; x) has a branch cut for x < 0; it is purely imaginary in the range —4 < x < 
and it is a complex number, with non vanishing real part, for x < —4. 



^The more extended set introduced in [TT] will be discussed in the following sections. 



The GHPLs with weight 2 or bigger are defined as repeated integrations of the weight 
functions in Eqs. (l2T]|22|) : 



G{a,w;x)= dtg{a;t)G{w;t), (26) 

Jo 

with the exception of G(0^; x), defined as: 

G{0^;x) = ^\og^{x). (27) 

Such a set of functions obeys (by construction) the shuffle algebra of Eq. ([7]) and all the 
important properties of the HPLs. 

The analytic properties of the functions defined in Eqs. ( 126)1271) derive from the properties 
of the logarithm and of the square root. 

2.2.1 Analytic Continuation 

When X > 0, every GHPL belonging to the set considered in the previous section is real 
and the only possible divergence is a logarithmic divergence in x = 0. 

For negative x, since the logarithm and the square root have a branch cut for negative 
argument, we must choose how to approach the cut. In order to do that, we give a vanishing 
imaginary part to the variable x. Let us choose the following prescription: 

X ^ -x' -iO-^ . (28) 

The region x < is divided in three sets, depending on the value of x'. The analytic 
continuation has to be done in each region differently. 

1. For —1 < X < 0, the imaginary parts come from log (x) and from the square root, 
that becomes purely imaginary: 

log(x) -)■ log(-x' - iO"*") = log(x') - ZTT, (29) 

1 1 ?■ . ^ 

^ .. , ., . .. =^^f^^- (30) 



v/x(x + 4) v^(-x' - i 0+)(-x' - i 0+ + 4) v^x'(4 - x') ' 

2. For —4 < X < —1, also the log (x + 1) gives an imaginary part: 

log (X + 1) ^ log (-X' - Z 0+ + 1) = log (x' - 1) - ZTT . (31) 

3. For — oo < X < —4, an additional imaginary part comes from the logarithm log (x + 4), 
while the square root becomes real again: 

log(x + 4) -^ log(-x'-iO+ + 4) = log(x'-4) -ivr, (32) 

i 1 , , 

^ .. , . ., . .. =-^^^^^- (33) 



v/x'(4 - x') v^(x' + iO+)(-x'-iO+ + 4) v^x'(x' - 4) ' 



2.2.2 Numerical evaluation 

The numerical evaluation of tlie GHPLs considered in the previous section can be done in 
principle using basically the same strategy of the HPLs. One focuses on x = 0, extracts the 
logarithmic behaviour using the shuffle algebra and then expands the remaining analytic 
functions (with no zeroes in the rightmost weight). Then, using suitable transformations, 
one relates the basic interval around x = to the rest of the real axis. 

However, the actual implementation of this strategy is quite cumbersome. Instead, it 
turns out to be convenient to transform from the beginning the GHPLs with square roots 
in the weights into a combination of GHPLs with linear weights using a set of variable 
transformations that will be discussed in the following section. The advantage of doing so 
lies in the fact that there exist fast and precise public numerical routines, that allow for the 
evaluation of generalized polylogarithms with generic linear weights [19]. The latter can be 
used to evaluate the GHPLs belonging to the set discussed in the last section, or, more in 
general, to the wider set introduced in [TI] . 



3 Linearization 

The presence, in an analytic result, of GHPLs with square roots together with linear weights, 
is due to the structure of the thresholds and pseudo-thresholds of the corresponding Feyn- 
man diagrams. 

Let us consider, for instance, the QED corrections to the vertex diagrams representing 
the decay of a photon into an electron-positron pair. The particle content reduces to 
massless photons and massive electrons/positrons. The threshold for the production of the 
electron-positron pair is at s = Am"^, while the pseudo-threshold lies at s = 0. Let us look at 
the differential equations with respect to s, for the solution of the corresponding Mis. The 
structure of the thresholds and pseudo-thresholds emerges in the homogeneous part with 
terms such as 1/s and l/(s — 4m^), that are also present in the non-homogeneous part (see 
for instance [22]). The solution of the homogeneous equation contains the inverse square 



root 1/ ^^/s{s — Am"^). The particular solution, then, comes from repeated integrations of 



1/a/s(s — Am?)^ 1/s, and l/(s — Aw?) terms. In the Euclidean region ijp' = — s > 0) the 
solution can be expressed in terms of GHPLs of the variable x = p^ /im? = —s/im?, with 
weights — r, —4, and 0. 

We can get rid of the square root (the weight — r) using the following quadratic trans- 
formation of variable: 

e ' ^ v/^T4 + yi' ^'^^ 



where ^ G C, |^| < 1, while x G M. In fact, we have: 

1 e 



(x + 4) (e + 1)'' 

1 ^ i 

^x{x + A) (e + l)(e-l) 



(35) 
(36) 



Moving from x to ^, the integration measure changes as follows: 

,,^/^ fa+l)(;-l) ,,,. (37) 

Ji T 

and every GHPLs reduces to a combination of repeated integrations of the simpler weight 
functions defined in Eqs. (II]). As a consequence of that, the set of weights {— r, —4,0} is 
transformed into the set { — 1,0, 1}, and the GHPLs are transformed into the usual HPLs 
defined in section 12.11 

Let us consider, now, a more complicated problem, in which zero- and multiple-mass cuts 
are present at the same time. This is, for instance, the case of the electroweak corrections to 
leptonic or hadronic processes in which the lepton and quark masses are neglected and only 
the vector boson masses are considered different from zero. In this case, the homogeneous 
part of the differential equations for the corresponding Mis contain terms as 1/s, l/(s — 
Am?), 1 / ^/s{s — 4m?) , together with terms as l/(s — m?). Therefore, the weights — r, —4, 
— 1, and are present at the same time. In this situation, it is more difficult to get rid of 
the square root. If we require, for instance, that the weights belong always to the set of real 
numbers, Wi G M, there is no transformation of variable that could linearize all the weights 
{— r, —4, —1, 0} at the same time. However, if we relax this constraint, we can move from 
the set with square roots and linear weights to a set of only linear weights, using the change 
of variable (134|) . 

Using Eq. ( 1M1) . the old weight functions are transformed into: 

with 

c= = e*3, c= = e*3, (42) 

where c and c are the two primitive sixth roots of the unity. Then, combining the integration 
measure, Eq. f l37|) . with the Eqs. f l38] - [4T|) . the original GHPLs with square root in the weight 
are transformed into 

G(-r, w;x) = [ dt g{-r; t) G(w; t) = - [ dr]- G(w; t(r/)) , (43) 

Jo Ji V 

^ ^ G{w;t{v)), (44) 



G(-4,w;x) = f dtg{-A-t)G{w-t)= f dr] (-- 
Jo Ji \ V 

G{-l,w;x) = / dtg{-l]t)G{w]t) = 
Jo 



v + i 



dv (-- + ^ + -^) G(w; t{r])) , (45) 

1 V^ rj — c f] — cj 

G'(0,w;x) = / dtg(Jd]t)G{yf;t)= [ drf { + ^— ) G(w;t(r/)). (46) 

Jo Ji \ V V-^J 

Therefore, the set {— r, —4, —1,0} of weights with square roots, has been transformed into a 
new set, with only hnear weights: {—1, 0, 1, c, c}. This new set, contains the original HPLs, 
discussed in section 12.11 and new GHPLs with complex weights c and c. The latter, have 
branch cuts in the complex x plain, starting at x = c, c respectively. At weight 1, they are: 

r dt 

G{c-x) = / = log (x - c) - log (-c) , (47) 

Jo t-c 

r dt 

G{c; x) = z = log {x -c) - log (-c) . (48) 

Jo t-c 

We can summarize the linearization procedure as follows: 

1. We base our analysis in the region in which x > (the formulas will be afterwards 
analytically continued in the region x < 0, if necessary). 

2. We transform the integration variable in the new variable t], on which we integrate 
from 1 to ^. Troubles with the integration in rj = 1 can occur, due to possible singular 
behaviours. Since such singularities can occur only from the weights in the variable 
X, we avoid the possible logarithmic divergence in rj = 1 using the shuffle algebra and 
extracting the trailing zeroes in x. The logarithms so found, G{On]x), are directly 
rewritten as l/n\ log" (x) and, then, straightforwardly transformed in the variable ^ 
using the relation log (x) = 2 log (1 — ^) — log (^). 

3. We linearize the GHPLs of weight 1. 

4. Weight-by-weight we proceed to the linearization of the GHPLs with weight 2 and 
bigger, integrating over the new integration measure the corresponding linearized 
GHPL times the corresponding linearized weight function. 

As an example, we give here the expressions of the linearized GHPLs with weight 1. 
G{0;x) can be converted directly in the new variable ^, sincqj: 

G(0; x) = log (x) = 2 log (1 - - log (0 = 2^(1; - ^(0; • (49) 

Using the relations in Eqs. ( 143] - [46|) . we have: 

G(-r;x) = -/ ^ = -G(0;O, (50) 

Ji V 

G{-4;x) = [ dv(-- + ^)=-2\og{2) + 2G{-l;0-G{0;0, (51) 

A V ^ v + ^J 
f^ f I 1 1 

G'(-l;x) = / dr] [ — + + 



1 \^ rj — c ^ — c^ 

G(c; 1) - Gic; 1) + G(c; + G(c; - G(0; • (52) 



Note again the different sign with respect to the weight +1 in the Remiddi-Vermaseren notation [5]. 
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It is worth to notice that the hnearization algorithm generates some constants, i. e. the 
hnearized GHPLs evaluated in ^ = 1. In many cases, these constants have a representation 
in terms of known transcendental constants. In general, however, this is not true. For the 
purpose of the numerical evaluation of the GHPLs with square roots in the weights using 
existing C++ routines, these constants can be left as they are. In fact, the routines provide 
a fast and accurate numerical evaluation in every point, and then also in ^ = 1. In our 
particular case we havqj: 

G(c;l) + G(c;l) = 0, (53) 

such that 

G(-l; x) = G(c; + G(c; - G(0; • (54) 

Knowing the expressions of the linearized GHPLs with weight 1 and linearized weight 
functions, we can proceed with the linearization of the GHPLs at weight 2. If we consider, 
for instance, the function G{0, —1; x) we have: 

G(0,-l;x) = / ^G(-l;t), 
Jo '^ 

= J^dv (-i + ^) [G{c; v) + G{c; v) - G(0; r^)] , 

= C(2) - G(0, c; - G(0, c; + ^G(0; 0' - 2G(1; 0^(0; 

+2G(0, 1; + 2G{1, c; + 2G{1, c; • (55) 

In the same way one can proceed for higher weights. Explicit formulas for the weight-2 and 
weight-3 GHPLs involved in the NLO electroweak corrections for the production of a Higgs 
boson in gluon fusion are provided in appendix |X] and appendix |Bl respectively. 

3.1 Analytic Continuation of the Linearized GHPLs and their 
Numerical Evaluation 

The analytic continuation of the linearized GHPLs is less complicated than the one concern- 
ing their original form. In fact, while the variable x ranges from oo to 0, the corresponding 
variable ^ is real and positive and it ranges from to 1. When x becomes negative, but in 
the range — 4 < x < 0, 

X ^ -x' - i 0+ , < x' < 4 , (56) 



^ becomes imaginary: 



y/x + A- a/x _ ^ a/4 - x' + iy/x' j2<t> 



e= ' V ^ C= , ' ^L— = p^^'P (57) 



where 



I X TT 

arctan \ , < < — . (5J 

— X 2 



^The constants (and also the GHPLs with complex weights in the actual expressions for the Mis) should 
appear always in a way such that their sum is real, as it has to be, since we are in the Euclidean region. 
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Finally, when x' ranges from 4 to oo, ^ becomes real again: 

C^r^^^-^. (59) 

and it ranges from 1 to 0. We must, therefore, discuss three regions. 

1. For < X < oo, we have < ^ < 1. The original GHPLs are real. The linearized 
GHPLs contain functions that are manifestly real, as the ones with weights —1,0, 1, 
but also functions that are complex: those that contain the weights c and c. However, 
the GHPLs containing the weights c and c appear in the formulas always in pairs (for 
instance G{—\,c\^ + G(— l,c;^)), in such a way that, although the single GHPLs 
of the pair are complex, their sum is real, since the imaginary parts are equal and 
opposite. The numerical evaluation of such GHPLs can be done straightforwardly 
using the VW routines presented in [19]. 

2. For —4 < X < 0, ^ is a pure phase, ^ = e*^''^. In order to evaluate numerically the 
GHPLs in this region, we have to notice that the VW routines, while allowing the 
use of complex weights, do not provide the possibility of evaluation of GHPLs with 
complex argument. However, for the GHPLs with non-trailing zeroes the following 
general formula holds: 

G{W\,W2-, ■■■■,VOn\X) = G{\W\,\W2-, ■■■■,y^'Wn\\x) , AgC, (60) 

as it will be discussed in the next section. Extracting the trailing zeroes and then 
choosing^ 

A = -, (61) 

X 

the GHPLs under consideration are transformed in GHPLs of real argument, ^ = \^ 
and complex weights, {±e~*^''^, 0, e~*^ sJ}^ that can be evaluated using again the 
VW routines. 

3. For — oo < X < —4, ^ is again real and we are back to the case explained in the first 
point. 

4 GHPLs of a Complex Variable 

In this section, we consider the case in which the GHPLs have to be evaluated in the 



complex plain. Therefore, x is complex from the beginning^. We are particularly interested 
in the following situation. Let us suppose that the dimensionless variable x is indeed a ratio 
between two physically meaningful variables: a squared momentum and a squared mass: 

2 

x = —^= ^ , 62 



® Actually, it is sufBcient to choose A — e""'"^^^). 

^The case in which x is real, but the corresponding reduced variable ^ is complex, was already discussed 
in the previous section. 
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with a/s the cm. energy of a certain process. This is, for instance, the case of the cor- 
rections presented in [21], but it is a quite general assumption. If the particle to which 
the mass m belongs is an unstable particle, its width F is going to play an active role in 
the determination of the corresponding physical observable. Consequently, the parameter 
X becomes complex, since we should now consider 



where 



"'" (m-^r/2)2~ M2''"^' ^^^^ 



M^ = m^ — , and = arctan < ^ J- . (64) 

4 rn? - 



In the non-physical region, in which s < 0, we have from Eq. fl63|) that 

7^e(x) > , and Xm(x) > . (65) 

The variable ^ defined in Eq. (1M|) . correspondingly, is also complex. The definition of the 
GHPLs does not change, except from the fact that now the integration is over a curve in 
the complex plain. Since the functions are analytic in the region defined by Eq. (1651) . the 
value of the GHPL does not depend on the path. In this region we have: 



_ Vx + 4 - y/j _ y/ri + xi + 4 - ^/r^^+xl + i {^/rl - Xi - 4 - y/ra - Xi) 
y/x + 4: + y/x A/ri +XI+4: + ^72 + xi + i (Vn - Xi - 4 + ^r2 - Xi) ' 



where: 



xi = 7^e(x) > , (67) 

X2 = Xm(x) > , (68) 



n = J(xi + 4)2 + xi, (69) 



r2 = \jx\^xl. (70) 

Therefore, we are in the situation in which we have to evaluate GHPLs with linear complex 
weights as functions of a complex variable ^. 

Let us consider a generic GHPL, Giwx^w^-, ...,Wn',x) in the case in which Wi,x G M. If 
no trailing zeroes are present, we can define a non- vanishing real parameter A G M, such 
that the following scale invariance holds: 

G{wi, W2, ..., Wn, x) = G{Xwi, Xw2, ..., Xwn] Ax) . (71) 

The demonstration of Eq. (I7ip can be done by induction. It is trivially verified for n = 1 
{X,Wi 7^ 0). In fact: 

/•Xx 

G{Xwi;Xx)= dtg{Xwi;t), (72) 

Jo 

and moving to the new integration variable r = t/X, we have: 

fX fX 

G{Xwi;Xx)= / Xdr g{Xwi; Xr) = / dr g{wi;r) = G{wi;x) . (73) 

Jo Jo 
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Let us suppose it is verified for n = i. For n = i + 1 we liave: 

G{Xwi-i-i, Xw; Xx) = / dt g{Xwi+i;t)G{Xw;t) = / Xdr g{Xwi; Xr)G{Xw; Xr) , 

Jo Jo 

PX 

= / drg{wi;r)G{w;r) = G{wi+i,w,x). (74) 

Let us suppose, now, tliat Wi,X,x G C For tlie weiglit 1 we have (remember tliat we 
are considering tlie case in wliicli |A|, \wi\ ^ 0): 

G{Xwi; Xx) = / dz g{Xwi; z) = / — , (75) 

Jo,7 Jo,7 z - Xwi 

where 7 is a path in the complex plain connecting the origin, 2; = 0, to the point Xx = 

\X\\x\e^^'^^^^'^''^"'^^^^'' = |A||a;|e**-'^^"^-'. If we rescale the integration variable by the real number 
|A||x|, we have 

GiXw,; Xx) = j^ ^ ^. _ ^^^i^+w,) ' C^^) 

where ^ = |w7i|/|x| and Wi = arg{wi). The path 7' connects the origin and the point on 
the circle of radius 1 with argument (A + X). Let us define 71 the path along the radius 
from the origin to e**^^"^"^) . F = 7' — 71 is a closed path that we suppose not to include the 
pole z' = ^e*^^"*"^^-*. The integral along the path F vanishes for the Cauchy's theorem. The 
integral over the radius can be rewritten as a one-dimensional integral of real variable with 
the substitution t = z' exp {—i{A + X)). Therefore: 

='''""^ dz' r dt 



G(Awi; Ax) = j^ ^^ ^; _ ^e.(A+m) = | t - ^e^^^-^) ^ ^^"^'Z'^' ^^ ' ^^^^ 



On the other hand, we have also: 



i(X) 



dz' r^ dt 



'^K5a:) = y^^^ -—-^^ = J^ ————- = Giwi/x-,l), (78) 

thus, 

G{Xwi;Xx) = G{wi;x). (79) 

Note that, in the end, for our purposes, we can just use Eq. (175]) . 

Let us suppose, now, that the rescaling is verified for n = i. For n = i + 1 we have: 

G'(Awi+i,Aw;Aa;) = / dz g{Xwi+i; z) G{Xw, z) = ^ _ (;(Aw; |A||x|z') , 

^0,71 J 0,^1^ s^ 

1 J. 

^ dt 

, ^ _ ^g^(H/,+,_x) G{w/x; t) = G{wi+i/x, w/x; 1) . (80) 
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Choosing A = 1 in Eq. fl5U|) . we can demonstrate that 

G{wi+i, w; x) = G{wi+i/x, w/x; 1) , (81) 

and, therefore 

G{\Wi+i, Aw; Ax) = G{wi^i, w; x) . (82) 

Using Eq. (18 ip . we can employ the numerical routines provided in [IS] for the evaluation 
of the GHPLs. In fact, now the GHPLs have complex weights (ratios of the original weights 
Wi and the variable x), but real variabl^J, equal to 1. 

Let us consider again our set {— r, — 4, — 1, 0}, and see what happens in the different 
regions. The analytic continuation from the non-physical s < region to the physical 
region in which p^ — )■ — s — iO~^, with s > 0, corresponds to the transformation x — )■ —x', 
where, now, x' G C and it is defined as follows: 

a;' = x;+^x^ = ^e^^ (83) 

Correspondingly, the variable ^ becomes C, -^ C^ with C G C defined as follows: 

^ yi^^-^/^^ ^ ^r[ -x[+A- y/r2 -x\-% (y/r^ + x^ - 4 - ^r^ + xQ 
V4 - x' + v^^a/ ^r\ - x^ + 4 + ^Jr^ -x\-i (^r^ + x'^ - 4 + ^r2 + Xi) ' 

where now 

r-'i = y^(4-xi)2 + xi. (85) 

Note that C^ does not have anymore modulus 1, as it was the case of real x shown in Eq. fH2|) . 

5 Generalizations and Additional Weights 

In this section, we enlarge the set of possible weights in order to cover the GHPLs needed for 
the analytic expressions of the Mis in [11]. The goal is to be able to describe the following 

set: 

{-\-r, -r, -4, -1, 0, 1, 4, r, 1 + r, c, c} , (86) 

where the additional weight functions (not introduced in the previous sections) are defined 
as follow^: 

^(4;x) = ^, (87) 

g{r-x) = / (88) 

a/x(x — 4j 

^(l + r;x) = , (89) 

y x(x — 4)(x — 1) 

^(-l-r;x) = , . (90) 

Vx(x + 4)(x + l) 



1° Note that it is sufficient to divide by e' ""'^i^') 

^-^Note the difference in sign in the definition of g(4; x), g{r; x), and ^(1 + r; x) with respect to [TT] . 
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The guidelines sketched in this section can be used for other, more comphcated, sets. 

It is first worth to notice that the possible weights listed in Eq. ( 186|) do not appear all 
together at the same time. The appearance of a particular weight in a GHPL depends on 
the cut structure of the relative Feynman diagram. In the Mis presented in [TI] we cannot 
have, for instance, the weights r and — r at the same time in the same GHPL. The same 
happens for the pair (c, c) with the square roots r or — r. Actually, the structure of the Mis 
in [H] is such that we are concerned effectively with three different subsets, that form each 
a closed base. They are: 

{-1, 0, 1, c, c} , {-1 - r, -r, -4, -1, 0} , {0, 1, 4, r, 1 + r} . (91) 

The three subsets do not mix with each other and they can be linearized (once and for all) 
using different variable transformations. 

5.1 The set { - l,0,l,c,c} 

The GHPLs belonging to this set can be evaluated straightforwardly with the help of the 
routines in [19] without any further variable transformation. In the case in which the 
variable x is complex, we just have to use the scale invariance of the GHPLs, as explained 
in section HI 

5.2 The set {- 1 -r,-r, -4, -1,0} 

This set contains the weights treated in section [3l {— r, —4, —1, 0}, with a small enlargement 
due to the weight (—1 — r). Note that this enlargement is totally painless, since the new 
weight (—1 — r) transforms in the same set of linearized weights { — 1, 0, 1, c, c}. In fact, 

r r dt 

G(-l-r,w;x) = / dtg{-l-r;t)G{w;t)= /-^-— , ^ G{w;t), 

Jo Jo Vt(t + 4)(t + l) 

= ^^ f dv (— —) G(w; tiv)) . (92) 

Note that the GHPL G(— 1 — r, w;x), which is real for x > 0, is written as a difference 
of the two complex GHPLs: (j(c, ...;^) and G(c, ...]■)]). This difference is indeed complex, 
since the two GHPLs have the same real part but opposite imaginary parts. The factorized 
i in Eq. (p2|) makes in such a way that the combination is real. 

5.3 The set {0, l,4,r, 1 + r} 

These positive weights cannot be linearized with the change of variable in Eq. f l34|) . Instead, 
we must use the change of variable that was used in [23] : 

,= (i±^, ^^ f-flj . (93) 

^ v^ + V a; - 4 
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When X is positive and ranges from cxd to 4, the corresponding variable u ranges between 
and 1. When < x < 4, u becomes imaginary. Giving to x a negative vanishing imaginary 
part (anticipating the prescription for the continuation to the Minkowski region), we have: 

e^2<^, (94) 



(95) 

(96) 
(97) 



\/x — \/x — 4 , a/x — ^/x — 4 - 

UJ — ) UJ — 


-zO+ ^/x + iVA- 


- X 


a/x + y/x -A y/x + \/x - 4- 

where 


-iO+ y^-iV4:- 


- X 


/4-x 








Finally, when x becomes negative. 






X —)■ -x' - iO^ , 


x' > 0, 




we have 


;- v^ 




„ Vx' + 4 





yjx' + 4 + yx' 

and uj" ranges between 1 and when x' ranges from to oo. 

Moving from x to w, the integration measure changes as follows: 

Jo J -I '/ 

Using eq. fl93|) the old weight functions are transformed into: 



1 ?7 

t^=(;^l) 



g{A-t) = -^ = ^-^, (101) 



^(r;t) = =-- -^(^ -, (102) 

""^ yt(^^ (r7 + l)(r/-l)' 

1 772 

g{l + r;t) = , = — -4- -, (103) 

^ v/t(^^(t-l) (r/ + l)(r/-l)(r/ + c)(r/ + c)' 

where the complex numbers c and c were defined in section [3l 

Combining Eq. ( 198|) with Eqs. ( I99l - I103p . we have the following transformation formulas 
for the definition of the GHPLs: 

G'(0,w;x) = / dtg{^-t)G{^-t)= f dr] f -^ + ^] G(w;t(r/)), (104) 

Jo J-i \ V V + ^J 

G{l,w-x) = r dtg{l-t)G{y^-t)= r dr, (-^- + ^^ + ^-\ ^(w; t(r/)) , (105) 
Jo J-i \ V V + c V + cJ 

G{4,w-x)=j dtg{4-t)G{w-t) = j dr, (-- + JL\ G{w;t{r])), (106) 
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G{r, w;x) = dt g(r; t) G(w; t) = - dr] - G(w; t{r])) , (107) 

Jo J-i V 

Gil + r,w-x) = r dtg{l + r-t)G{w-t) = i^rdr^ ("^ _ ^"j G(w; t(r/)). (108) 
Jo 3 J_i \r]+c r]+cj 

The integration in t] deserves a further discussion. As in the case aheady presented in 
section [3l the point t] = —1 can be source of a non integrable singularity. However, the 
possible divergence in r^ = — 1 is connected to the original point x = 0, and then, ultimately, 
to the right-most weight in the GHPLs of x. It is sufficient, therefore, to extract the right- 
most trailing zeroes in x before the change of variable (I93p is applied, using the shuffle 
algebra. The functions G(0„; x) = l/n\ log" (x) can be directly transformed in the new 
variable u using the relation log (x) = 2 log (a; + 1) — log(ci;). The GHPLs that do not 
contain trailing zeroes in the right-most weights are regular in r] = —1 after the variable 
transformation. 

5.4 Mixed Weights 

Although the weights belonging to the different sets described above do not mix in the 
expressions of the Mis of [11], we can further extend the analysis and try variable trans- 
formations that linearize wider sets of weights. This can be done provided that we do not 
mix the square roots with different signs. For instance, it can be shown that the weights 
belonging to the set { — 1 — r, — r, —4, —1, 0, 1, 4} can be linearized at the same time, using 
the variable transformation in Eq. flM|) . Analogously, the set {—4, —1, 0, 1, 4, r, 1 + r} can 
be linearized with the help of the change of variable of Eq. fl93l) . 

6 Two-loop Light-Fermion contributions to the Higgs 
Production in Gluon Fusion 

In this section, we revisit the calculation of the NLO light-fermion electroweak corrections 
to the Higgs boson production in gluon fusion. 

In [21] these corrections were evaluated analytically, and the results were expressed 
in terms of GHPLs with square root in the weights. The numerical evaluation was done 
using real W and Z masses and with FORTRAN routines written ad hoa^\ The electroweak 
corrections appear to be very peaked at rriH ~ ^miy and rriH ~ ^mz because of the opening 
of the two corresponding thresholds. In this section we recompute the NLO-EW corrections 
using the VW routines employing complex values for the W and Z masses. As a result, the 
finite W and Z widths smear the peaks at the thresholds and resize the relative importance 
of the corrections in the region mn ~ 2mw, 2mz- 

Neglecting QCD corrections, the partonic production cross section, up to 2-loop level. 



^^ In [24], the remaining electroweak corrections due to the top quark were calculated as a Taylor 
expansion in mj^ / (Arn^y) . Finally, in |25j a numerical calculation with complex W and Z masses was done 
for the complete set of NLO electroweak corrections. 



has the following form: 



^i99^H) = ^^:^\g'^ + agf '^ 



512^2 



7r 



'EW\ ' 



(109) 



where Gp is the Fermi constant, as the strong coupling constant, and a the fine structure 
constant. 

The lowest order, Q^\ is due to one-loop diagrams with heavy quarks running in the 
loop. The dominant contribution comes from a loop of top, while the contribution of a 
b-quark loop is of the order of some percents of the previous one. The analytic expression 
of g^^ is: 



g- = J2 



q=t,b ^1 



2- 1 + 



4 



Xr 



G{-r,-r]Xg) 



(110) 



where Xg = —rn^j/m?, rriq is the heavy-quark mass (top or bottom mass), and the analytic 
continuation has to be taken considering a positive vanishing ttlh imaginary part: Xg — )■ 
—x' — iO+, where x' = m^j/m?. G{—r,—r;Xg) can be immediately transformed into a 
square logarithm of the variable ^, defined in Eq. (1341) . as for instance in Eq. (IllSp . 

The two- loop electroweak light-fermion contributions, gfl, to ^^V can be expressed as 



m-. 



72/ 



[mw — i^wf^Y 



yif 

where s^ = sin^ 9w, c^ = ^ — s'^, 



2-715^ m'jj 



2/5 7 
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^ 7 - o ^ + TT^ M^z) + 4:A,{xw) 



(111) 



771 



Xw 



H 
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{mw — iVw/'^Y ' 



xz 



H 



{mz - iTz/2Y 



(112) 



and 

AAx) 



-4 + 2(1 + -) ^(-1; x) + -^(0, -1; x) + 2 ( 1 + - ) ^(0, 0, -1; x) 
xj X \ xj 

+ ( 1 + ^ j [2^(0, -r, -r; x) - 3G'(-r, -r, -1; x)] - ^^(x + 4)|^G'(-r; x) 



+ 



x + 2 



x^ 



2G(-r,-r,-r;x) + 2G(-r,0,-l;x)-3G(-4,-r,-l;x) \. (113) 



The GHPLs with square root in the weights involved in Eq. (I113P are listed in section |3] 
and in appendix iBi 
Writing 

where ctq is: 



a{gg ^H) = ao{l + 6, 



If) ' 



ctq 



GpOis |^1Z|2 



512y2 7r 



gY, g'' = gl' + gl'. 



we have for 5, 



If- 



2a 



^if 



\Q 



\l\2 



\Re{g^') ne{gf}) + Xm(6;") Xm(^j 



(114) 
(115) 

(116) 
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Figure 1: 5// as defined in Eq. flll6p . 
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Table 1: 5// as a function of tlie Higgs boson mass imu in GeV). 



In Fig. [T] we plot bij computed with the VW routines and some numerical values are collected 
in Table [H The set of parameters used is the following: 



mi = 173.1 GeV, mft = 4.6GeV, mty^ 
mz = 91. 1876 GeV, F^ = 2.4952 GeV, 
a = 1/128, Gir = 1.16637 -10-^ GeV"^ 



80.398 GeV, Fvf = 2.141 GeV, 

0.23149 (117) 



sin^ 6'vK 



7 Conclusions 

In this paper we analyzed the set of GHPLs of a single variable containing square roots in 
the weights. After recalling the definition and basic properties of the HPLs, we introduced 
the GHPLs with weights belonging to the set {— 1 — r, — r, —4, — 1, 0, 1, 4, r, l-f r, c, c}. This 
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specific set of GHPLs appears in tlie analytic expressions of the Mis that enter into the 
calculation of the electroweak form factor [201 E] • 

One of the main observations of the paper lies in the fact that, once the weights are 
allowed to be complex, the GHPLs with square roots in the weights can be "linearized", 
i.e. expressed as a combination of GHPLs with linear weights. These linearized GHPLs 
are functions of a transformed variable, that is not unique, but can be properly chosen 
depending on the nature of the weights. The set {—1 — r, — r, —4, —1, 0, 1, 4, r, 1 + r, c, c} 
can be linearized, once and for all, with just two variable transformations. 

The other observation concerns the possibility of a fast and precise numerical evalua- 
tion of the linearized GHPLs using already existing numerical routines. In particular, the 
C++/GiNaC routines by Vollinga and Weinzierl [12] offer a well suited tool for this goal. 

Finally, the strategy for the numerical evaluation of GHPL presented in the paper 
is applied to the known case of electroweak light-fermion NLO corrections to the Higgs 
production in gluon fusion. We evaluate the GHPLs with square roots using the VW 
numerical routines. As a further refinement, while in [21] the corrections were evaluated 
neglecting the effects of the W and Z widths, we consider here the case of complex mw 
and mz, getting a more realistic result. It is worth to notice that the GHPLs with square 
roots allow for a very compact analytic expression of the results, which would be extremely 
lengthy if expressed in terms of the linearized GHPLs. 
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A Some Examples at Weight 2 

As simple examples, in this appendix we apply the procedure outlined in the paper to the 
following GHPLs at weight 2: G{—r, —r; x) and G{—r, —1; x). We find: 

G(-r, -r; x) = f J^ G{-r- 1) = - f '^ [-G(0; r/)] , 
Jo yjt{t + 4) Ji V 

= Ig{0;O\ (118) 

G(-r,-l;x) = r J^M -l- 1) = - t^ [G{c- v) + Gic; v) - G{0; v)] , 
Jo ^/t{t + 4) Ji V 

= -ic(2)-G(0,c;O-G(0,c;O + iG(0;O'. (119) 
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B Some Examples at Weight 3 



In this appendix, we provide the expressions for the 5 GHPLs at weight 3 containing square 
roots in the weights, involved in the corrections of [21]. We have: 



G(— r, — r, — r; x) 



dt 



G{0, — r, — r; x) 



^/t{tT^ 

-^G'(0;O^ 

"" df 
^G(-r,-r;t) 

^ 



G(-r,-r;t) 



dr] 1 
^2 



G{0;vY, 



1 \ V Tl-lJ 2 



G'(-r,0,-l;x) 



2C(3) - -G(0; if + G(0; O'^ll; + 2G(0, 0, 1; 
-2G(0;OG(0,l;O, 

G(0,-l;t), 



Vt(t + 4) 

1 V 
-2G{1- ri)G{^- 7]) + 2^(0, 1; r/) + 2^(1, c; r^) + 2^(1, c; r/) 



C(2)-G(0,c;r7)-G(0,c;r7) + -G(0;r/)2 



(j'(— r, — r, —1; x) 



-^C(3) + 27^1 - C(2)G(0; - ^^(0; 0' + ^(0, 0, c; 

6 

+G(0, 0, c; e) - 2G(0, 1, c; - 2G(0, 1, q 
+2G(0;OG'(0,l;O-2G(0,0,l;O, 

^i=G(-r,-l;t), 
^t{t + 4) 

1 



1 ^ 



C(2) + G'(0,c;r7) + G(0,qr/)-G'(0,0;r/) 



G(-4,-r,-l;x) 



^C(3) + k(2)G(0; - Jg(0; 0' + G(0, 0, q 
3 3 

+G(0,0,qO, 

t + 4 ^ ' ' ^' 
5 



(120) 



(121) 



(122) 



(123) 



-C(2) + G(0, q r/) + G(0, q r^) - G(0, 0; r^) 



df] [ 

1 V^ ^ + 1 

-|C(3) - 2i^2 - |C(2)G(-1; - 2G(-1; 0^(0, q 
D 3 

-2G(-1; 0^(0, q + G(-l; 0^(0; 0' + 2G(0, q -1; 
+2G(0, q -1; + ^C(2)G(0; - ^G(0; 0' + 2G(0, -1, q 
+2G(0, -1, q - 2G(0; 0^(0, -1; + ^(0, 0, q 
+G(0,0,qO + 2G(0,0,-l;O. 



(124) 
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In the formulas above, we introduced the two constants Ki and i^2- They have a 
cumbersome expression in terms of known transcendental constants, that we omit here. 
Their numerical value is known with infinite precision and it is: 

Ki = 0.278425076639727748441973590814.., (125) 

K2 = -0.152226248227607546589100778278... (126) 
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